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ABSTRACT 

We present a new mechanism for Type la supernova explosions in massive white dwarfs. The 
proposed scenario follows from relaxing the assumption of symmetry in the model and involves a 
detonation created in an unconfined environment. The explosion begins with an essentially central 
ignition of stellar material initiating a deflagration. This deflagration results in the formation of a 
buoyantly-driven bubble of hot material that reaches the stellar surface at supersonic speeds. The 
bubble breakout forms a strong pressure wave that laterally accelerates fuel-rich outer stellar layers. 
This material, confined by gravity to the white dwarf, races along the stellar surface and is focused 
at the location opposite to the point of the bubble breakout. These streams of nuclear fuel carry 
enough mass and energy to trigger a detonation just above the stellar surface. The flow conditions 
at that moment support a detonation that will incinerate the white dwarf and result in an energetic 
explosion. The stellar expansion following the deflagration redistributes stellar mass in a way that 
ensures production of intermediate mass and iron group elements consistent with observations. The 
eject a will have a strongly layered structure with a mild amount of asymmetry following from the early 
deflagration phase. This asymmetry, combined with the amount of stellar expansion determined by 
details of the evolution (principally the energetics of deflagration, timing of detonation, and structure 
of the progenitor), can be expected to create a family of mildly diverse Type la supernova explosions. 
Subject headings: hydrodynamics — instabilities - stars: interior — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae are one class of luminous stellar 
explosions. These are the predominant explosive events 
in old stellar environments such as elliptical galaxies. 
The ejecta of these objects are rich in intermediate mass 
and iron group elements. Explaining the nature of these 
objects is therefore critical for understanding galactic 
chemical evolution (jTruran fc Cameron! I1971D . These 
supernovae also are the key component of one method 
used to determine the histo ry of the Universe and probe 
the origin of dar k ener gy ([ Sandage fc Tammannl I199& 
jPerlmutter et, aTl Il993 Tfonrv et all 1200.1 iKnon et, all 

The work presented in this Letter builds on many pre- 
vious observational and theoretical contributions to the 
field of Type la SNe, and on advances in fluid dynam- 
ics, nuclear physics, and computational science. Cur- 
rent ideas about the Type la supernova explosion mech- 
anism follow from the original work of Arnett, Nomoto, 
and Khokhlov (IArnetflll96fl iNomoto. Sugimoto. fc Ned 
119761 lKhokhlovlll991h . who pioneered deflagrating and 
detonating massive wh ite dwarfs originally proposed by 
iHovle fc Fowlerl ()1960D as the core component of Type la 
supernovae. 

Despite decades of effort, these events remain 
an unsolved mystery. Proposed exp losion scenar- 
ios include w hite dwarf detonations (jArnettl Il969t 
iNomotol 11982ft . coalescing pairs of white dwarfs 
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(jWebbinkl 11984 llben fc Tutokovl 119841 . deflagrations 
or delayed detonations of mas sive whi te dwarfs 
(INomoto. Thielemann, fc Yokol [HH |Khokhlov| 
I1991D. and collapse in a strong gravitational field 
({Wilson fc Mathews] 120041) . None of these scenarios 
accounts for all the observed features of Type la super- 
novae. Some models produce energetic explosions but 
fail to explain the observed ejecta compositions, while 
others successfully produce the observed chemically 
stratified ejecta but require including ad hoc physics. 

Here we report the results of multi-dimensional hydro- 
dynamical simulations of the long-term evolution of a 
massive white dwarf following an essentially central igni- 
tion. The initial evolution results in a deflagration front 
and formation of a Ray leigh- Taylor unstable buoyancy- 
driven bubble that is accelerated to supersonic speeds 
on its way to the stellar surface. We follow the evolution 
beyond bubble breakout and observe formation of grav- 
itationally confined flow across the surface of the star. 
The flow is focused at the point opposite the breakout 
location, where the matter is compressed and heated, ig- 
niting a detonation wave just above the stellar surface. 
We find that the star expands substantially during the 
evolution. This expansion produces a density distribu- 
tion that, when overrun by the detonation wave, can be 
expected to result in ejecta that are strongly layered and 
rich in intermediate mass elements, as is typical of Type 
la supernovae. The modest asymmetry in the expan- 
sion of the star can also be expected to produce a mild 
amount of global asymmetry in the explosion. 

2. NUMERICAL MODEL 

The simulations presented here were performed with 
the ad aptive mesh refin ement hydrodynamics code 
FLASH ((Frvxell et al.ll2000l) . The numerical scheme in- 
cludes self-gravity solved using a multipole expansion. 
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We have established that to properly account for asym- 
metries in the mass distribution and to ensure momen- 
tum conservation in the model, at least 3 multipole mo- 
ments are needed in the expansion. In the simulations 
reported here, we used 10 multipole moments. The rest 
of the physics modules are identic al to those used in our 
previous study ([Calder et all2()(M) , in particular the evo- 
lution of the deflagration front is computed with a flame 
capturing scheme with energy release accounting for ca r- 
bon, magnesium, and silicon burning ( K hokhlovll200lt) . 

The computational domain is a two-dimensional re- 
gion in cylindrical geometry covering a region from 
— 16,384 km to 16,384 km in the z-direction and ex- 
tending up to 16, 384 km in radius. We employed an 
outflow-only boundary condition at all domain bound- 
aries, except along the symmetry axis at r = where 
we used a reflecting boundary condition. The adaptive 
mesh was used to track flow features at a maximum ef- 
fective resolution of 8 km (corresponding to an effective 
grid resolution of 2048 x 4096). The simulation was ex- 
ecuted at Courant number of 0.6, and the time step was 
limited to a maximum of 4 x 10 -4 s. 

The initial conditions consist of a 1.36 solar mass 
isothermal white dwarf composed of equal amounts of 
carbon and oxygen with a temperature of 3 x 10 7 K. The 
progenitor model was mapped to th e computational gri d 
following the procedure employed in lCalder et al.l(|2004l ). 
The nuclear flame was initiated as a small spherical re- 
gion of completely burned material placed in hydrostatic 
equilibrium with its surroundings. The ignition region 
centered at (r, z) = (0, 12.5) km had a radius of 50 km. 

3. RESULTS 

Figure [I] depicts several moments in the evolution of 
the density of the star from the point of bubble breakout 
to just prior to the ignition of the detonation. The rapid 
ascent of the bubble toward the stellar surface acceler- 
ates the stellar material located just above the bubble. 
This piston-like behavior results in the formation of a 
bulge filled with high-pressure high-momentum nuclear 
ash. During breakout, bubble material is expelled mostly 
radially, while the high pressure of the burned material 
accelerates the surface layers laterally (Fig. ^a)). This 
material races along the stellar surface, followed by the 
magnesium-rich bubble material. Both remain gravita- 
tionally confined to a relatively thick w 1000 km layer 
(Fig. IHb)). At t « 1.8 s, this flood of surface material 
converges at the point opposite to the bubble breakout 
location, forming a conical compressed region bounded 
by the shock. This structure stretches down the symme- 
try axis beginning at z « — 3 x 10 8 cm (Fig. EJc)). 

At this point in the evolution, conditions in the shocked 
region approach the detonation regime. By t = 1.85 s, 
material upstream of the confluence region has density 
of ~ 10 4 g cm -3 and moves at velocity 9,500 km s _1 . 
Downstream of the shock, the density of the nuclear 
fuel is ~ 5 x 10 4 g cm -3 and the temperature reaches 
^ 1.4 x 10 9 K. Due to a mild density gradient present 
in the flooding material (Fig. Efa)), the post-shock den- 
sity slowly increases with time while the temperature re- 
mains relatively constant. At t « 1.93 s, the post-shock 
conditions are suitable for igniting the nuclear fuel: the 
density exceeds 1.7 x 10 6 g cm -3 and the temperature 
is « 2.2 x 10 9 K. The detonation point can be seen as a 



slightly over-pressured region located near the symmetry 
axis at (r,z) « (0,-3.35 x 10 8 ) cm (Fig.^b)). 

The detonation wave born in the region above the stel- 
lar radius will sweep through the white dwarf, which un- 
derwent a substantial evolution from the moment of the 
ignition of the deflagration. During its ascent, the rising 
bubble displaced about 5% of the stellar mass. This mass 
displacement, combined with the pressure wave caused 
by nuclear energy release, leads to expansion of the star. 
Thermal expansion is present from the onset of the defla- 
gration with that part of the star nearest the deflagration 
experiencing relatively stronger thermal expansion. The 
global evolution of the stellar matter is, however, pri- 
marily in response to the softening of the gravitational 
potential caused by mass displacement. This expansion 
has a mostly radial character and becomes effective only 
after bubble breakout. 

Figure |3| shows the evolution of density during the 
course of the simulation. Shown is the amount of mass 
in several density intervals. The material with density 
> 1 x 10 8 g cm -3 will undergo burning to iron peak ele- 
ments, while material with density below 3 x 10 T g cm -3 
will burn into intermediate mass elements. Several points 
can be noted here. There is almost no stellar expansion 
prior to the bubble breakout (t < 0.9 s). The break- 
out is followed by relatively fast expansion of the densest 
material (p > 1 x 10 9 g cm -3 ), which is succeeded by 
uniform expansion (identified by a simultaneous increase 
in the expansion rates at t ~ 1.2 s). This process contin- 
ues until the ignition of the detonation, with the most 
rapid expansion at densities between 1 x 10 7 g cm -3 
and 1 x 10 8 g cm -3 . At the time of detonation the 
amount of mass at the densities characteristic for produc- 
tion of intermediate mass elements, p < 3 x 10 T g cm -3 , 
is 0.35 Mq. The amount of mass at densities above 
1 x 10 8 g cm -3 , required for production of the iron peak 
elements, is 0.71 Mq. 

4. DISCUSSION AND CONCLUSIONS 

We presented a gravitationally confined detonation 
(GCD) mechanism for Type la supernova explosions. 
The basic components of the proposed scenario are a 
rising deflagrating bubble expelling a small amount of 
stellar matter, the associated stellar expansion caused 
by the shallower potential well, the flood of stellar mate- 
rial across the surface following the bubble breakout, and 
a detonation in an gravitationally confined environment. 

In the GCD scenario, not imposing any asymmetries 
in the initial conditions is of paramount importance. A 
deflagration front is born very close to the white dwarf 
center. Such an essentially central ignition is more prob- 
able than the idealized conditions adopted in standard 
deflagration or delayed detonation models, primarily be- 
cause the central region of the star is convective. This 
type of ignition results in a rising deflagrating bubble ac- 
celerated by buoyancy to supersonic speeds on its way 
to the stellar surface. The transonic phase of the bub- 
ble's rise is accompanied by the formation of a bow shock 
ahead of the bubble that compresses and heats the nu- 
clear fuel. Our attempts to associate this region with a 
possible transition to detonation failed. We discovered, 
however, that the flood of the expelled surface layers that 
follows the bubble breakout remains confined to low radii 
above the star. This flood races around the star and is 




Fig. 1. — Evolution of the white dwarf following an essentially central ignition. Density in log-scale is shown for the innermost 
approximately one-sixth of the simulation domain, (a) The bubble breakout (t = 0.9 s). The material is expelled radially, and the high 
pressure of the burned bubble material produces a lateral acceleration of the outer layers of the star, (b) The bubble material expands 
above the star as the fuel-rich surface layers reach the equator in their race across the stellar surface (t = 1.4 s). Note that most of the 
surface layers are closely confined to the star, (c) Focusing of the fuel-rich streams just prior to the detonation [t = 1.9 s). Note the 
presence of the dense conical region stretching down along the symmetry axis beginning at (r, z) = (0, —3 x 10 8 ) cm. 




Fig. 2. — Evolution of the surface flood following an essentially central ignition of the white dwarf. The lower half of the star is shown in 
the figures, (a) Density in log-scale at t = 1.85 s. Note the density gradient present inside the stream of fast-moving material flowing across 
the stellar surface. Also note that the interior of the star remains highly radially symmetric, and that only relatively low-density regions 
are perturbed by the surface flow, (b) Pressure in log-scale at the time of the ignition of the detonation (t = 2.005 s). The ignition point 
can be seen as a small over-pressured region located near the symmetry axis at (r, z) « (0, —3.25 x 10 8 ) cm. The shock wave preceding the 
bubble material can also be seen in the low-density stellar layers near (r, z) « (1.1 x 10 8 , —2.4 x 10 8 ) cm. The contour marks the position 
of the advancing front of burned bubble material. Vectors show the velocity field; the magnitudes of these velocities can be determined by 
comparing the length of these vectors with the length of the fiducial vector, which corresponds to 1 x 10 9 cm s _1 . 



ultimately focused into a hot, compressed, high density 
region located just above the stellar surface. Conditions 
in this region satisfy the criteria necessary for a detona- 
tion. 

One important aspect of the GCD mechanism is that 
stellar expansion is a natural consequence of the essen- 
tially central ignition of a deflagration. The flame re- 
leases energy and displaces mass, softening the gravita- 



tional potential well leading to expansion of the star. 
The expansion will slow down on a time scale compa- 
rable to the sound crossing time of the white dwarf as 
the star tries to reach a new equilibrium. Therefore, we 
expect that at still later times, if not for the fact that 
the detonation will completely disrupt the white dwarf, 
the initial expansion would be followed by contraction of 
stellar material and the star would oscillate. 
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Fig. 3. — Evolution of the mass distribution in the model white 
dwarf. Amount of material in solar masses inside select density in- 
tervals is shown as a function of time. Significant stellar expansion 
takes place only after bubble breakout (about 0.9 s after the defla- 
gration began near the center of the star). The expansion velocity 
increases linearly from the stellar center and leads to a steady de- 
crease in the density of the star at all radii. At the moment of 
the ignition of the detonation, almost half of the stellar mass has 
densities below 1 x 10 8 g cm -3 . 



Because of this pre-expansion, the detonation front 
born above the stellar surface will encounter den- 
sities similar to those found in models where pre- 
expansion res ults from a centrally ignited large-scale 
deflagration |Reinecke. Hillebrandt. fc Niemeverl 120021: 
iGamezo et al.ll2003ft . The estimate of nucleosynthetic 
yield for intermediate mass (iron peak) elements is a 
lower (upper) limit in view of the fact that the stellar 
expansion continues after the moment of the detonation. 
Determining the final yields, however, requires simulat- 
ing the detonation. The actual conditions across the det- 
onation wave, particularly the amount of compression, 
will influence the results. Also, the yields can be mod- 
ified by delaying or by accelerating the ignition of the 
detonation to create a family of Type la supernovae with 
diverse spectral characteristics. The exact timing of the 
detonation depends on several factors. The structure of 
the progenitor influences energy release by the deflagra- 
tion. It also affects the strength and mass of the stellar 
layers being pushed around the star. The radius of the 
progenitor regulates the time required by the streaming 
matter to reach the confluence point. All these factors 
determine the time available for stellar pre-expansion and 
will be a source of diversity in GCD models. 

Some properties of the proposed model are, how- 
ever, largely independent of the precise details of 
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the ignition of the detonation or stellar progenitor. 
The explosion will be powerful. All the stellar fuel 
will be consumed by the detonation. Despite the 
perturbation introduced by the deflagration, the star 
will retain most of its radially symmetric stratifica- 
tion. Therefore, the subsequent explosion will display 
characteristics well-known from one-dimensi on al in- 
vestigations (iNomoto. Thielemann. &; Yokoil 119841 
iHoflich. Wheeler, Thielemann! I1998D . In particular , 
we expect the distribution of nucleosynthetic products 
in velocity space to agree with the observed layered 
structure of Type la supernova ejecta. 

The model also naturally admits certain asymmetries. 
The deflagration consumes only about 5% of the stellar 
mass. We expect a similar level of variation in the result- 
ing spectra and luminosities, in agree ment with degr ee of 
diversity present in the observations ([Li et al.ll2001l) . In 
addition, because the stellar shape is distorted by the 
rising bubble and the formation of the detonation on one 
side of the star, we expect a noticeable asymmetry in the 
stellar ejecta. These orientation effects might be respon- 
sible for peculiar events such as SN 1991T ([Filippenkol 
I1997D . Because the gross properties of observed Type 
la supernovae can be accounted for in the GCD model, 
detailed observations will be required to verify the pro- 
posed mechanism. One possibility is to obtain informa- 
tion about the degree of asymmetry in these explosions, 
using precision polarimetric measurements. 

In this short communication, we have not presented 
simulations that support all of our predictions. Rather, 
we have presented a logical sequence of events that nat- 
urally leads to a Type la supernova explosion. Several 
of these steps need to be carefully studied. In particular, 
the proposed detonation mechanism bears many similari- 
ties to the process of confined fusion studied in terrestrial 
laboratory experiments, which is notoriously difficult and 
prone to instabilities. For this reason, the early stages of 
the detonation should be studied very carefully. The ex- 
pected computational demands and challenges are severe 
and clearly approach limits of feasibility. Despite the fact 
that such a study lies in the future, we are confident in 
the basic components of the proposed GCD mechanism. 



The authors thank Alexei Khokhlov and Peter Hoflich 
for their comments. Contributions from Natasha 
Vladimirova, Ed Brown, Jim Truran, and the FLASH 
Code Group made this work possible. This work is sup- 
ported in part by the U.S. Department of Energy under 
Grant No. B341495 to the Center for Astrophysical Ther- 
monuclear Flashes at the University of Chicago. 



Khokhlov, A. M. 2001, ApJ, submitted ( astro-ph/0008463 ) 
Knop, R. A., et al. 2003, ApJ, 598, 102 
Li, W., et al. 2001, ApJ, 546, 734 
Nomoto, K. 1982, ApJ, 253, 798 

Nomoto, K., Sugimoto, D., & Neo, S. 1976, Ap&SS, 39, L37 
Nomoto, K., Thielemann, F.-K., & Yokoi, K. 1984, ApJ, 286, 644 
Perlmutter, S., et al. 1999, ApJ, 517, 565 

Reinecke, M., Hillebrandt, W., & Niemeyer, J. C. 2002, A&A, 391, 
1167 

Sandage, A., & Tammann, G. A. 1993, ApJ, 415, 1 



Gravitationally Confined Detonation 



5 



Tonry, J. L., et al. 2003, ApJ, 594, 1 Wilson, J. R., & Mathews, G. J. 2004, ApJ, in press. 

Truran, J. W., & Cameron, A. G. W. 1971, Ap&SS, 14, 179 
Webbink, R. F. 1984, ApJ, 277, 355 



